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^ ABSTRACT: The KATRIN experiment aims at the direct model-independent determination of the 

i-C average electron neutrino mass via the measurement of the endpoint region of the tritium beta 

decay spectrum. The electron spectrometer of the MAC-E filter type is used, requiring very high 
^ stability of the electric filtering potential. This work proves the feasibility of implanted 83 Rb / Kr 

calibration electron sources which will be utilised in the additional monitor spectrometer sharing 
the high voltage with the main spectrometer of KATRIN. The source employs conversion electrons 
of 83m Kr which is continuously generated by 83 Rb. The K-32 conversion line (kinetic energy of 
17.8keV, natural line width of 2.7 eV) is shown to fulfill the KATRIN requirement of the relative 
energy stability of ±1.6ppm/month. The sources will serve as a standard tool for continuous 
monitoring of KATRIN's energy scale stability with sub-ppm precision. They may also be used 
in other applications where the precise conversion lines can be separated from the low energy 
. ^ spectrum caused by the electron inelastic scattering in the substrate. 
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1. Introduction 

The KArlsruhe TRItium Neutrino (KATRIN) experiment represents a next-generation tritium j3 
decay experiment designed to perform a high precision direct measurement of the electron neutrino 
mass m(v e ) with the sensitivity of 0.2eV/c 2 (90 % C.L.). In the /3 decay experiment the observable 
is the weighted squared mass flU] 

3 

"* 2 (Ve) = £|£/e;| 2 m 2 (v 7 -), (1.1) 
7=1 

where U denotes the unitary mixing matrix and Vi, V2, V3 are the mass eigenstates. The square root 
of eq. fll.ip is often called the average electron neutrino mass. KATRIN is a successor experiment 
of the neutrino mass experiments carried out in Mainz (Germany) [Q] and Troitsk (Russia) [01 
which set the upper limit on m(v e ) of about 2eV/c 2 [@L ED- Therefore, the aim of the KATRIN 
experiment [U, 0, 0] represents the improvement of the neutrino mass sensitivity by one order of 
magnitude. 

For the observation of a non-zero neutrino mass signature in sub-eV range in the endpoint 
region of tritium j3 spectrum (decay energy Q = 18.6 keV) the methods of high resolution electron 
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spectroscopy are necessary together with a high luminosity and very low background level. For this 
purpose the upcoming KATRIN experiment [0] uses a windowless gaseous tritium source and two 
successive electrostatic retardation filters with the magnetic adiabatic collimation (so-called MAC- 
E filters) [0, 0D- However, the three aforementioned requirements are not the only ones which are 
connected to the challenging realisation of the KATRIN experiment. A possible instability of the 
energy scale of the main KATRIN spectrometer is one of the main systematic effects: The princi- 
ple of the MAC-E filter technique relies on the knowledge of the electric retarding potential which 
is experienced by the electrons on their path through the spectrometer. Owing to the designed 
sensitivity, the precise knowledge of the retarding potential at every stage of the measurement is 
inevitable in the KATRIN experiment. Besides the use of the state-of-the-art equipment for a direct 
measurement of the high voltage (HV), several very stable calibration electron sources, based on 
atomic/nuclear standards, will be utilised in KATRIN. One of the electron sources will be con- 
tinuously measured by an additional MAC-E filter spectrometer (so-called monitor spectrometer 
which is an upgraded spectrometer of the former Mainz Neutrino Mass Experiment) to which the 
HV will be applied, corresponding at the same time to the filtering potential of KATRIN. This way 
a redundant twofold monitoring system will be formed, increasing the confidence in the stability of 
the main spectrometer energy scale. In this work the feasibility of solid electron source based on 
the metastable isotope 83m Kr is tested. In this type of source the process of internal conversion of 
83m Kr is utilised, where 83m Kr is continuously generated by 83 Rb implanted into the substrate. The 
monitoring task of KATRIN demands the relative energy stability of 8E /E < ±1.6 ppm / month of 
the K-32 conversion line (kinetic energy of E = 17.8 keV, natural line width of T = 2.7 eV [Q]). 

This article is organised as follows: After briefly introducing the KATRIN experiment, in sec- 
tion we explain the necessity of the energy scale stability and its monitoring. The requirements for 
an electron source for such monitoring are discussed and it is shown that the implanted 83 Rb/ 83m Kr 
source may fulfill these requirements. In sections |3] and |] the test measurements carried out at the 
MAC-E filter spectrometer in the Institut fur Physik, Universitat Mainz, are described. Firstly the 
experimental set-up and evaluation methods are introduced (section 0) and later the experimental 
results obtained with the first four samples of the implanted 83 Rb/ 83m Kr source are summarised 
(section |j). Finally, in section |5] we discuss our results and give an outlook on the further develop- 
ment of the calibration electron sources for KATRIN. 

2. Natural calibration electron standard for KATRIN 
2.1 The MAC-E filter technique 

A sensitive search for the neutrino mass by measuring the endpoint region of the tritium /3 decay 
spectrum requires an electron spectrometer with high luminosity as well as high energy resolution. 
Both these basic requirements are fulfilled in the concept of magnetic adiabatic collimation with 
an electrostatic filter, abbreviated as MAC-E filter. This principle, firstly introduced in the 1970s 
[0] and early 1980s flTO, ED, Q2L was later adopted for the use in the neutrino mass experiments at 
Mainz (Germany) [@] and Troitsk (Russia) fl3fl. 

Basically, a spectrometer of the MAC-E filter type is a vacuum vessel containing a cascading 
system of cylindrical HV electrodes and two superconducting solenoids which are placed on both 
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ends of the vessel. The solenoids generate a highly inhomogeneous magnetic field guiding the 
electrons along the magnetic field lines from their origin in the source to the detector. The magnetic 
field is symmetrical with respect to the central plane (so-called analysing plane) of the spectrometer. 
The minimum field strength B m [ n at the central plane is reduced by several orders of magnitude with 
respect to the maximum B ma x occurring at the centre of the solenoids (also called pinch magnets). 
The magnetic gradient force transforms most of the transverse energy E± of the electron into the 
longitudinal motion. This is represented as a smooth collimation of the momentum vector with 
respect to the magnetic field direction. As far as the magnetic field B varies relatively slowly along 
one cyclotron loop of the electron the energy transformation takes place adiabatically, i.e. the 
magnetic moment pL keeps constant 1 [QJ, 

E± 

jU = — = const. (2.1) 

Thus, the j8 electrons, isotropically emitted from the source of a given area, are transformed into 
a broad beam of electrons flying almost parallel to magnetic field lines. This parallel beam of 
electrons is running against an electrostatic potential formed by the set of cylindrical electrodes. All 
electrons with enough energy to pass the electrostatic barrier are reaccelerated and collimated onto 
a detector, all others are reflected. Therefore the spectrometer acts as an integrating high-energy 
pass filter. Varying the electrostatic retarding potential allows the measurement of an electron 
spectrum in an integrating mode. 

The energy resolution of a MAC-E filter follows from eq. fl2.1|). In the extreme case where 
the total kinetic energy of the electron at the starting point, E stan , is given in the form of transverse 
energy, the adiabatic transformation according to eq. fl2.i| ) will result in a small remainder of trans- 
verse energy left over at the analysing plane. This maximum amount of non-analysable energy 
defines the theoretical resolution AE, 

B ■ 

AE = (E±) max = E stait - — — . (2.2) 

^>max 

Therefore, the resolution is basically only limited by the minimal ratio B m [ n /B max of the magnetic 
fields (at the central plane and at the centre of the solenoids) that can be realised experimentally. 
In order to suppress electrons which have a very long path within the tritium source and therefore 
exhibit a high scattering probability, the source is placed in a magnetic field B s < B max . Due to the 
magnetic mirror effect the maximum accepted starting angle d™£ of the electrons is restricted to 

s ^ = arcsin,/^. (2.3) 

V #max 

For an isotropically emitting electron source the transmission function of an ideal MAC-E 
filter can be expressed analytically [0, 0] as a function of the retarding potential U < and the 
magnetic field strengths B s , B max and B m i n . 

1 Since the electrons generated in tritium j5 decay reach a maximum value of the relativistic parameter y = 1 .04 « 1, 
the non-relativistic approximation may be used for this brief explanation of the principle. 



-3- 



2.2 Motivation for continuous monitoring of the KATRIN energy scale 

The complex set-up of the KATRIN experiment, spanning over 70 m in length and currently be- 
ing built and commissioned at the Karlsruher Institut fur Technologie (KIT), can be described as 
follows: The high luminosity windowless gaseous tritium source ensures a constant count rate of 
j3 electrons which are magnetically guided to the tandem of MAC-E filter spectrometers. The so- 
called pre-spectrometer rejects the large low energy part of the tritium j8 spectrum which does not 
carry any significant information for the determination of m(v e ). Thus, only the electrons origi- 
nating from the last « 300 eV below the tritium j8 spectrum endpoint of Eq = 18.6keV enter the 
succeeding so-called main spectrometer. Here the high precision energy analysis takes place. The 
electrons which pass the analysing potential are counted in the detector (silicon PIN diode array of 
148 pixels with a total area of 63 cm 2 ) with the energy resolution of ~ 1.5 keV for 18.6keV elec- 
trons. The kinetic energy of the /3 electrons will be analysed on the energy scale which is defined 
by the difference between the tritium source potential and the analysing retarding HV potential of 
— 18.6kV of the main spectrometer. 

The fluctuations of the energy scale represent one of the six main systematic uncertainties of 
the KATRIN experiment. The KATRIN project requires that none of the individual systematic 
effects gives the uncertainty contribution of more than about Am 2 yst ■ = 0.007 eV 2 /c 4 , thus the total 
systematic uncertainty (where further minor contributions are included as well) with respect to the 
observable m 2 (v e ) is estimated as A/n 2 yst = 0.017eV 2 /c 4 {J. The source potential will be a few 
hundreds Volts which can be measured precisely with the help of commercial devices. On the other 
side, the precise measurement of the HV is generally less reliable as instabilities can occur. The 
relationship between the instabilities of the retarding HV and the systematic uncertainty of m (v e ) 
can be derived by approximating the j8 spectrum by a simplified form in its endpoint vicinity and 
considering the HV fluctuations, introduced to the measurement. Assuming the HV fluctuations as 
a Gaussian distribution with a given width a one arrives at UJ, IT3H 

Am^(v e )c 4 = -2<r 2 . (2.4) 

This relationship shows that an unrecognised fluctuation of the retarding potential in the form of 
a Gaussian broadening of the width a imprints on the /3 spectrum in the manner of reducing the 
neutrino mass squared, i.e. the introduced HV fluctuation affects the /3 spectrum as a shift towards 



negative m 2 (v e ) values. Finally, comparing eq. (2.4) with the KATRIN requirement, the upper limit 



on the unrecognised fluctuation of the HV scale is obtained: 

Am 2 €V (v e ) < 0.007 eV 2 /c 4 => a < 0.059eV. (2.5) 

The upper limit on a stated in eq. ( [231 ) can be interpreted as the requirement on the HV scale 
stability of 

A(HV) abs = ±60mV (2.6) 
which represents the relative stability of 

A(HV) rel = ±3.2ppm (2.7) 

of the retarding potential of 18.6kV corresponding to the tritium jS spectrum endpoint Eq. The 
indicated HV stability must be fulfilled over each KATRIN experimental run spanning about two 
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months 2 . A separate value of m 2 (v e ) will be determined from each run and the resulting neutrino 
mass will be calculated as an average value. Evidently, successful stabilisation of the HV scale 
over more than one run would allow to directly combine the individual data series and to examine 
the j6 spectrum at higher statistics. The total live time for achieving the aforementioned neutrino 
mass sensitivity amounts to three years Hjfl. 

For the purpose of an ultra-precise monitoring of the energy scale stability in KATRIN a 
twofold monitoring system is currently under construction in the KIT. Firstly, the HV of — 18.6 kV 
applied to the main spectrometer will be measured directly using a specially developed high preci- 
sion HV divider together with a state-of-the-art digital voltmeter. Secondly, a suitable calibration 
source of monoenergetic electrons will be continuously measured in the third MAC-E filter spec- 
trometer of KATRIN, the aforementioned monitor spectrometer. This spectrometer will share the 
HV with the main spectrometer. Any change of the measured energy of the monoenergetic electrons 
can indicate a problem in the direct measurement of the HV. This way an important redundancy 
can be achieved by relying not only on high precision devices, but also on the natural standard of 
atomic and nuclear physics. 

The direct measurement of HV with the sub-ppm precision will be accomplished with the help 
of two ultra-stable custom-made high precision HV dividers flM . Oil developed for KATRIN in 
the Institut fiir Kernphysik, Universitat Mimster, in cooperation with the Physikalisch-Technische 
Bundesanstalt (PTB) Braunschweig. The HV divider scales the HV of —18.6 kV down to low volt- 
age of about — 10 V which can be measured by a commercial digital voltmeter with high precision. 
The voltmeter is regularly calibrated with the help of a commercial 10 V DC reference which, in 
turn, is calibrated against PTB's losephson voltage standard IITHl . The HV divider concept is sim- 
ilar to the one used before for producing the PTB's standard divider MT100 flT7ll . although very 
different types of resistors were used in these two concepts. Both KATRIN dividers K35 (designed 
for HV up to 35 kV) and K65 (max. HV of 65 kV) were extensively tested and calibrated against 
the PTB's standard HV divider MT100 and also against each other. The fact that two such devices 
are available for KATRIN is of a great importance with respect to their mutual cross-checks and 
overall redundancy. The dividers K35 and K65 were also recently utilised [ITEIl as the reference for 
calibration of the HV installation at the ISOLDE facility at CERN. To our knowledge the dividers 
K35 and K65, together with the PTB's divider MT100, represent the three worldwide best dividers 
in the given HV range. 

2.3 Requirements for an electron standard 

On the contrary to e.g. gamma spectroscopy, where various gamma ray standards fITPIl are com- 
mercially available, the situation in the electron spectroscopy is much less favourable. To our 
knowledge there exists no commercial standard of a solid electron source with the energy stability 
required by KATRIN. The consideration of the HV stability constraint in KATRIN expressed in 
eq. ( [2.5D leads to the following requirements for a suitable calibration electron source intended for 
the application at the monitor spectrometer: 

2 Tritium run times longer than two months will be avoided due to the time needed for the regeneration of special 
cryogenic components in the experimental set-up. 
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1 . The source should provide a well defined discrete electron line, ideally with the energy close 
to the endpoint of the tritium j3 spectrum. 

2. The electron line energy should be stable enough to check the relative stability of the main 
spectrometer retarding HV set by eq. ( [2.7D over the time period of at least two months. It 
should be noted that the foreseen application of the source is the recognition of the HV 
scale fluctuations on the time scale of minutes to months. On the other hand, the HV scale 
fluctuations of the order of seconds cannot be recognised with this method. Such short-term 
fluctuations can be assessed with the help of sophisticated HV ripple probes. 

3. It may seem unnecessary to state such a stability requirement with respect to an electron 
calibration standard as e.g. the gamma ray standards naturally provide stable gamma lines. 
Similarly, this holds also for an electron source in gaseous phase. However, in the case 
of the conversion electron emitted from a solid source the binding energy of the electron 
is influenced by solid state effects which directly translates into the kinetic energy of the 
ejected electron. Therefore, the solid source should provide stable chemical environment for 
the emitted electrons. 

4. The electron line should be sharp, i.e. the intrinsic energy spread should be as small as 
possible. In the case of a natural standard this requirement translates into a very small natural 
line width T of the electron line. 

5. In case of a natural standard the half-life of the source should be suitable for its use during 
the KATRIN experimental run of two months. Ideally, the half-life should span more than 
one run in order to allow for combined analysis of several data taking runs with tritium. 

6. For the calibration purposes only the zero-energy-loss electrons are suitable, in other words 
the electron which undergoes an inelastic collision in the source is lost from the useful elastic 
peak and therefore not usable for the calibration. The count rate of the zero-energy-loss elec- 
trons should be high enough in order to provide good statistics in a reasonable measurement 
time. 

2.4 Implanted 83 Rb/ 83m Kr source 

An electron source based on the internal conversion of 83m Kr was already used in the past for 
the calibration purposes in Mainz Neutrino Mass Experiment H2DH and other neutrino experiments 
flZTl U3\ . In KATRIN the 83m Kr source will be applied in three different physical states [fj: as a 
gas, as a condensed film of the sub-monolayer thickness and as a solid source. Regardless of the 
physical state of the given source, the parent isotope 83 Rb with the half-life of 7\/ 2 = 86.2(1) d 
[H3] will be exploited as 83 Rb conveniently serves as a generator (via electron capture) of the short- 
lived isomer 83m Kr with Ti/ 2 = 1.83 h (1T9T1 . The decay scheme of 83 Rb is depicted in figure [jj. The 
conversion electrons produced in electromagnetic transitions from the krypton nuclear levels above 
the isomeric state are not usable for our purpose due to their high energies. However, the isomeric 
state decays via a cascade of suitable low energy transitions of 32.2 and 9.4 keV which possess 
high intensity of conversion electrons. The first transition (32.2 keV) with the multipolarity E3 is 
especially highly converted: The total internal conversion coefficient (ICC) amounts to a tot = 2010 
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Figure 1. The decay scheme of 83 Rb to 83 Kr (left scheme) and a zoom into the decay of metastable 83m Kr 



to Kr (right scheme). The Rb decays by a pure electron capture with the branching ratio of 77.9 % into 
the short-lived isomeric state 83m Kr. Detail on the right side shows its 32.2keV (intensity of 0.0358(45) %, 
multipolarity E3) and 9.4 keV (5.86(134) %, M1+E2) gamma transitions. The schemes are taken from iPl l. 



IE3I1 . The second transition (9.4 keV) is practically a pure Ml transition with a small E2 admixture 
(mixing parameter 8 = 0.0130(8)) and the total ICC amounts to a tot = 17 flZBTl . 

The gamma ray energy of the 9.4 and 32.2 keV transitions was determined by gamma ray spec- 
troscopy with semiconductor detectors as E 32 = 32 151.7(5) eV [El and £9.4 = 9405.8(3) eV []25D, 
respectively. The most interesting conversion line of 83m Kr regarding the application in KATRIN 
is the one created by the internal conversion of the 32 keV transition on the K shell as its kinetic 
energy is only by 0.8 keV lower than the tritium endpoint. The precision of the electron binding 
energy E^' ga& (K) = 14327.26(4) eV of this core shell for gaseous krypton was increased in the 
reevaluation \UM . The kinetic energy of an electron emitted from atomic shell i is roughly given by 

JMO^r-JC""®, (2-8) 

where E y is the gamma ray energy and E^' gas (i) is the electron binding energy (related to the 
vacuum level) of the shell i of a free atom. 

The solid type of 83m Kr source was introduced in the framework of the KATRIN experiment 
as a convenient alternative to the proven principle of the condensed 83m Kr source (CKrS) [0, EZD 
which is based on earlier work of the Mainz Neutrino Mass Experiment flZQl . The concept of such 
a solid source is based on the fact that the 83 Rb generator of 83m Kr lies directly in the source itself. 
Thus, the count rate of the source is driven by the half-life of 83 Rb. The systematic measure- 
ments of the energy stability of the 83m Kr conversion lines, carried out at the Mainz MAC-E filter 
spectrometer, were started with 83 Rb/ 83m Kr sources prepared by vacuum evaporation of 83 Rb onto 
carbon substrate. In the course of these measurements, which are reported elsewhere EEL EJ E3h 
it was realised that the sources fulfil the requirements for the time stability of the conversion line 
energy but exhibit only moderate retention of 83m Kr of about 15 %. The retention of 83m Kr repre- 
sents the portion of 83m Kr atoms retained in the solid matrix of the source and thus useful for the 
calibration purposes. In addition, the vacuum-evaporated sources were shown to be very suscepti- 
ble to surrounding conditions (humidity, vacuum). The idea to implant the 83 Rb ions into a suitable 
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substrates was based on the early work IBTH where the 57 Co ions were implanted at the energy of 
500 eV into cleaved surfaces of natural graphite crystals. 

The implanted 83 Rb/ 83m Kr source should be a very promising alternative to the CKrS and 
even superior in terms of vacuum requirements and ease of use. However, certain change of the 
electron binding energy is to be expected due to solid state effects 3 . This can be illustrated in 
accordance with the approach used in [fJll . The kinetic energy of the conversion electron emitted 
from the atomic shell i of the 83m Kr atom placed in a metal host and measured by the spectrometer 
reads H 

■^kirf 1 (0 = By + ^rec, 7 ~ £rec, e (0 _ E^' ™ Pl (0 ~ (j>spec ~ ^source) , (2-9) 

where imp (i) denotes the electron binding energy related to the vacuum level. The term E TeC} r 
denotes the energy of the recoil atom after gamma ray emission and is equal to 0.007 eV and 
0.002 eV for the 32.2 keV and 9.4 keV transitions of 83m Kr, respectively. The energy E rec , e (i) of 
the recoil atom after the emission of the conversion electron from shell i is smaller than 0.22 eV and 
0. 12 eV for the two mentioned transitions in 83m Kr, respectively. The work functions sp ec, ^source of 
the spectrometer electrode and source, respectively, are important for relating the measured kinetic 
energy to the retarding potential of the spectrometer. The quantity impl (/) can be expressed in 
the form 

£ b T imP '(0 = E™< g3S (0 - AE^(i) , (2.10) 

where AE^(i) > represents the solid state correction. This correction related to the Fermi level 
of the metal host reads 

AEgr(i) =A£S(0 + Source- (2-11) 
In this case the eq. ( |2T9"1 ) can be rewritten into 

EZ Vl (})=Er+Er e o,y-E Tec>e ({)- (^ p, (0-AEgr , (i)) (2-12) 

The samples investigated in this work were prepared via ion implantation at the ISOLDE 
facility at CERN. The 83 Rb ions were implanted into polycrystalline foils of gold and platinum 
of high purity. A solid uranium carbide target was used in combination with tantalum surface ion 
source, as rubidium, being an alkali metal, is easily ionised. In order to maximise the portion 
of zero-energy-loss electrons a low implantation energy of 30 keV was chosen, corresponding to 
the minimal implantation energy available at ISOLDE at that time. In one case, the ion beam 
was decelerated by retarding HV down to 15keV. In all the cases the target foils were kept at 
room temperature and perpendicular to the ion beam direction. The implantation took place under 
relatively poor vacuum conditions of « 10~ 4 mbar. Prior to the implantation, the foils of 12 mm 
diameter were cleaned in reagent-grade methanol by ultrasonics. In addition, the beam focusing 
was optimised using a non-radioactive Rb beam of high intensity hitting a paper replacing the foil. 
This was important especially for the implantation at 15keV, in which case the foil was placed 
about 1 m behind the usual place of beam incident. The option of sweeping the ion beam over the 

3 In literature the phenomenon of core level binding energy shifts of rare gases implanted in noble metals was studied 
both experimentally El El El EH and theoretically EE El El El- Nonradioactive rare gases (Ne, Ar, Kr and Xe) 
were implanted into polycrystalline foils of high purity and the XPS method was utilised for measuring the absolute 
electron binding energies. These were then compared with the binding energies measured for free atoms. 
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foil was not utilised during our collections. After the completion of collections the foils were stored 
on air in special holders with plastic caps preventing any wipe-off. It should be emphasised that no 
annealing treatment of any sample was carried out, thus the damage of the lattice by the irradiation 
was not healed. However, as the collections took place at room temperature it can be assumed 
that the lattice underwent self-annealing to some extent. It is expected that after its implantation, 
83 Rb behaves as Rb + and bonds to the environment of the Au or Pt lattice. Some kind of diffusion 
cannot be excluded, but due to its reactivity 83 Rb will probably stay on its original place. Some 
strong bonding of 83m Kr in the solid phase is not expected as within « 10~ 10 s after the decay of 
83 Rb the electron shells are reorganised. However, as the metastable 83m Kr exists for rj 2 h before 
the conversion electron is emitted, the Au or Pt lattice surrounding the 83m Kr atom might influence 
the conversion electrons via weak bonds. 

The first 83 Rb collection resulted in the 3.3 MBq source denoted as Pt-30 (30keV 83 Rb ions 
hitting the Pt target). The polycrystalline foil of thickness of 40 used for this collection was 
not a standard commercial foil. Its purity was determined as 99.7 % with the help of the X-ray 
fluorescence method, the 0.3 % impurity being identified as rhodium. For the second series of 
collections the commercially available foils 4 of thickness of 25 ptra, purity of 99.99+ % and temper 
quality "as rolled" were used, with the impurities specified at the level of tens of ppm or less. As 
the first implanted source Pt-30 exhibited promising results concerning the energy stability of the 
83m Kr conversion lines, the same implantation energy and foil element was used once more (but the 
foil was purer in the second case). This led to the sample designated as Pt-30-2 (4.9 MBq of 83 Rb). 
However, the two Pt-30 sources were not strictly identical due to the different levels of impurities 
and implanted doses. In order to investigate other possibilities of the 83 Rb implantation the sources 
Au-30 (3.3 MBq of 83 Rb in gold foil) and Pt-15 (1.9 MBq of 83 Rb, implantation energy of 15 keV) 
were produced as well. The amount of 83 Rb present in the sources right after the implantation 
was deduced from measurements of the 83 Rb activity by gamma spectroscopy carried out about 
5 months later. 

2.5 Distribution of the ^Rb atoms in the implanted sources 

The areal distribution of 83 Rb atoms over the surface of the four sources was studied iWTTl with the 
help of a Timepix detector P3H . The scans of the 83 Rb activity are shown in figure ||. It was found 
that the distribution of 83 Rb in the plane of the three sources implanted with a 30keV beam can be 
roughly represented by a two-dimensional Gaussian distribution with a given widths a,b (FHWM). 
On the other hand, the source Pt-15, produced with the help of the retarding electrode, exhibits 
a circle-segment-like shape. From the knowledge of the implanted 83 Rb activity and the activity 
spot area the mean implanted dose can be calculated. In all the cases the dose laid in the range of 
« (1-6) x 10 14 ions/cm 2 (see table [[]) at a rate below 2.5 x 10 9 ions/s (corresponding to a current 
of 400 pA). From the activity distributions (see figure 0) it can be estimated that the maximal areal 
dose power did not exceed 220 pt W / cm 2 so that any heat load caused by the beam can be neglected. 
A dose of the order of 10 14 ions/cm 2 is generally considered as a limit above which the host lattice 
damage takes place [B3H . As expected, in the case of the well focused sources Pt-30-2 and Au-30 
the implanted dose was the highest. The retention of 83m Kr in the implanted sources was found to 

4 The Au foils AU000341/49 and the Pt foils PT000240/55 from the company Goodfellow E3 were used. 
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Pt-30 Pt-30-2 Pt-15 Au-30 

Figure 2. The areal distribution of the 83 Rb activity in the implanted 83 Rb/ 83m Kr sources. The red circles 
mark the source holders of the inner diameter of 1 1 mm. The figure is based on IH-UI . 



be > 90 % [FTTH . It may seem surprising that the retention is not 100 %, however, it was pointed out 
in API PT5H that a diffusion of the implants to the surface might occur. 

The distribution of 83 Rb atoms along the depth of the foils was assessed with the help of the 
SRIM code [ffoTl . Simulations were carried out with the following parameters: The implantation 
energy was set to 15keV (source Pt-15) and 30keV (all other sources) and the incident angle 
(relative to the source foil normal) was set to zero. There is no temperature parameter in the 
SRIM simulation 5 . Polycrystalline structure of the foils is the default and only option in the SRIM 
code describing well our situation. Due to poor vacuum conditions surrounding the foils during 
the implantation, there was certainly a considerable amount of water and hydrocarbons adsorbed 
on the foil surface. This fact was reflected in the simulations by introducing an additional layer 
above the metal foil. For the sake of simplicity, this "rest gas layer" was set to be a pure carbon 
layer of the thickness of 3 nm which was assumed to be a good approximation of the surrounding 
conditions. Moreover, the sputtering effect, produced by 83 Rb ions impinging on the surface, is 
negligible as the ions possess too high energy and penetrate too deep into the metal foil. This 
statement is supported by the simulation results: The sputtering yield of 30 keV Rb ions was found 
to be in the range of 10-20 atoms (Pt or Au) per Rb ion which, together with rather low 83 Rb ion 
beam intensities of « 2 x 10 9 ions/s at ISOLDE, means that less than one monolayer of the rest 
gas layer was sputtered during the implantation. 

The implantation ranges, obtained after the simulation of 1.5 x 10 5 of ions in each case, are 
depicted in figure |[ It can be seen that the portion of 5-10 % of the incident 83 Rb ions gets stopped 
in the rest gas layer which represents a non-metal environment. As the foils were not cleaned by ion 
sputtering or any other means after the collections, a certain portion of 83 Rb is assumed to remain 
in this layer. The simulated values of the mean projected range R p and the longitudinal straggle 
AR P (standard deviation of the range distribution) are also summarised in table |I[ Knowing the 
values Q and AR P , via the approximation [B7I| 

Q = V2n-C p -AR P (2.13) 

one can estimate the peak concentration C p of 83 Rb. In all four cases the peak concentration lies 
in the range of (1-6) x 10 20 ions/cm 3 . Considering the density of Pt or Au atoms as 6.62 x 

5 The self-annealing effects, occurring at room temperature, are not included in the SRIM calculation 1 14611 . 
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Table 1. Overview of the implanted 83 Rb/ 83m Kr sources. A stands for the 83 Rb activity right after the source 
production. The values a,b are the widths (FWHM) of an elliptic shape of the activity spot. Knowing the 
values A,a,b the mean implanted dose Q was calculated. The simulations carried out with the SRIM code 
delivered values of the mean projected range R p and its standard deviation AR p . The peak 83 Rb concentration 
C p is estimated on the basis of Q and AR p . Considering the density of lattice atoms, the relative peak atomic 
concentration c p of 83 Rb is also calculated. In the last column the measured retention of 83m Kr is listed. 
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Figure 3. The range profiles of Rb ions implanted into gold and platinum, simulated with the help of the 
SRIM code: Pt-30 (green points), Pt-15 (magenta) and Au-30 (blue). The dashed line at the depth of 3 nm 
denotes the boundary between the metal foil and the carbon layer mimicking the rest gas layer which was 
included in the simulation. The maxima of the distributions were individually scaled to unity. 

10 22 /cm 3 or 5.90 x 10 22 /cm 3 , respectively, the estimated 83 Rb peak atomic concentrations are 
found in the range of 0.2-0.8 %. 



3. High resolution conversion electron spectroscopy with the MAC-E filter 
3.1 Experimental set-up 

The samples of the implanted 83 Rb / 83m Kr sources were investigated with the electron spectrometer 
of the MAC-E filter type in the Institut fur Physik, Universitat Mainz. This MAC-E filter, shown 
in figure |], was previously used as the j8 electron spectrometer in the former Mainz Neutrino Mass 
Experiment [@J] until 2001 and later it was upgraded to the high energy resolution of 0.9 eV at 



During the Mainz Neutrino Mass Experiment the spectrometer was typically operated at the energy resolution of 
4.8eVat 18.6keV. 




Figure 4. The photograph of Mainz MAC-E filter set-up in 2009: (i) solenoid A, (ii) solenoid B, (iii) spec- 
trometer vessel wrapped in thermal insulation, (iv) auxiliary coils for adjusting the spectrometer resolution, 
(v) auxiliary perpendicular coils for magnetic flux positioning and Earth's magnetic field compensation. The 
electron sources and the detector were placed in the bores of the solenoids B and A, respectively. 

17.8 keV. 

The filter consisted of an ultrahigh vacuum tank (length of ~ 4 m, diameter of ~ 1 m), massive 
conical and cylindrical inner electrodes (stainless steel, DIN 1.4429), two "massless" wire elec- 
trodes and a magnet system. The vacuum tank and the conical massive electrodes on opposite sides 
were grounded. The wire electrodes were conductively connected with the massive electrode and 
a negative analysing voltage of up to —35 kV was applied to them. The work function spec of the 
central analysing electrode was measured as spec = 4.4(2) eV [0]. The magnet system comprised 
two superconducting solenoids A and B and two sets of auxiliary air coils. The first one, consisting 
of four coils, was used for adjusting the field fi m i n , i.e. the energy resolution of the spectrome- 
ter. The second set of two windings, perpendicular to each other, allowed to create a magnetic 
field in any transversal direction. This was used for compensating for misalignments of the source 
and detector relative to the symmetry axis (z direction) of the spectrometer. Working positions of 
the source and detector were on the z axis inside the bores of the corresponding superconducting 
solenoids. The magnetic field in the centre of solenoids amounted to 6.014 T. The distance between 
the centres of the solenoids was 4.02 m. The minimum magnetic field was set to B m ; n = 0.3 mT 
in order to achieve the spectrometer energy resolution of AE = 17.8keV • 0.3mT/6T = 0.9eV at 
17.8keV, see eq. flZJ). 

The electron sources were placed in a holder, see figure |3[ The source holder was mounted on 
a long rod fixed at the flange. The flange was attached to the x-y-z table which allowed to centre any 
of the four sources onto the spectrometer axis. The x-y-z table was connected via a large diameter 
bellows to the spectrometer so that the positions of the sources could be changed without breaking 
the vacuum. The movement in the z direction was motor-driven with the precision of 1 mm. The 
source could be positioned in the range of 12.5-70 cm relative to the centre of the solenoid B. The 
movement in the x and y directions was done manually with the precision of 0. 1 mm. The elevation 
amounted to ±1 cm for both directions. The samples were electrically isolated among each other 
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Figure 5. The photograph of the source holder allowing to investigate up to four samples without breaking 
the vacuum. On the right side the detail of the four foils (each in its own sample holder) with the implanted 
83 Rb/ 83m Kr sources is shown. The gold foil is recognisable in the top right position, the other three foils are 
the platinum foils. The sample holder in the bottom right position was made from aluminium, whereas the 
other holders were produced from gold plated copper. The inner and outer diameters of each sample holder 
were 11 and 14 mm, respectively. 

and also relative to the spectrometer. Each sample was connected to the vacuum feedthrough by a 
wire which allowed to ground or bias the sample, the maximal applicable voltage was 3 kV. This 
way it was possible to separate each electron source energetically (see section The usual 
measurement position of the source was chosen as z = 21 cm from the centre of the solenoid B, 
corresponding to B s = 1.75 T and 6™n = 32.7°, cp. eq. (|Z.3|). Larger acceptance angles 0™ a ™ were 
avoided in order to limit the energy losses of electrons starting inside the source foil under large 
angles. 

The windowless Si PIN photodiode Hamamatsu S3590-06 with the effective depletion layer of 

0. 5 mm thickness [ffEU was used as the electron detector. Its active area of 9 x 9 mm 2 was mounted 
beneath a circular mask of diameter 9 mm. The detector set-up was mounted on a movable x-y-z 
table which allowed to adjust the detector inside the bore of the solenoid A with respect to the 
spectrometer z axis. The usual position of the detector was ~ 23 cm from the solenoid A centre, 

1. e. in the magnetic field of 1.06 T. The detector was placed on a LN2 cooled finger together with 
the preamplifier electronics. The energy resolution was 3 keV (FWHM) at 17.8 keV, the noise edge 
of the photodiode was found to be 6.4 keV. Thus, the measurements of the low energy conversion 
lines at 7.5 keV were still realisable but with losses in efficiency. The detector signals were ampli- 
fied by a custom-made charge sensitive preamplifier, shaped by a spectroscopy amplifier (Silena 
7614, shaping time 0.5 jits) and digitised by a 13-Bit CAMAC ADC (EG&G ORTEC AD413). 
Parallel to the detector signals the signals from a pulse generator (EG&G ORTEC 448, frequency 
100 Hz) were fed to the charge sensitive preamplifier in order to determine the dead time T of the 
system. The frequency stability of the pulser was regularly checked in such a way that the pulser 
signal was delivered directly to the spectroscopy amplifier instead of the preamplifier signal. In 
addition, a histogramming memory (EG&G ORTEC HM413, FERA-bus interface) was utilised, 
which significantly decreased the dead time. The dead time was determined to be z ~ 40 jUs in 
the event-by-event readout mode and T < 10 /is in the histogramming readout mode. The precise 
timing of the measurement was based on the /i-scaler (SPEC. 003, frequency 100 kHz) which was 
used together with the precision frequency standard (AR 750747, frequency 1 MHz). 

The high precision HV power supply (FuG Elektronik HCP- 18-350000, stability 2ppm per 
8h, reproducibility lOppm [fl"9H , voltage ripple 20 mV peak-to-peak, ppm stability warmup time 
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2 h [1501 ) delivered the HV to the massive and wire electrodes inside the spectrometer. The HV was 
also fed to the aforementioned high precision HV divider K65, scaling the HV down to 20 V range. 
The key property of a HV divider is the long-term stability of its dividing ratio M(t). In general, 
the long-term dependency of M(t) can be described as exponential [15011 which — after sufficient 
time has passed from the divider production date — can be simplified to a linear function 



where the long-term drift m is typically of the order of 0.1 ppm/month for a high precision device. 
The long-term stability of the dividing ratio M = 1818.1096(36) of the HV divider K65 was 
extensively studied in the Institut fiir Kernphysik, Universitat Miinster, over the time period of 
13 months [fl"5L BTTl . The dividing ratio drift was determined as m = —0.13(10) ppm/month, i.e. 
practically compatible with a zero drift. 

The low voltage in the 20 V range was measured by a commercial digital voltmeter (Fluke 
8508A, 20 V DC range, 7 ! /2-digits resolution mode [52TD . The voltmeter was regularly (once per 
1-2 days) calibrated with the help of the 10 V DC voltage reference (Fluke 732A, stability specified 
as 3ppm/year 152T1 but actually determined as < 0.3ppm/year H5T11 7 ) which was calibrated at the 
PTB. This way the offset and the drift of the scale factor (gain) of the voltmeter were monitored 
and taken into account. 

The experience gained at the PTB showed that the main uncertainty in the low voltage mea- 
surement comes from cable effects like the resistance and capacity of the coaxial cables used. In 
addition, the thermoelectric effect gives rise to thermo-voltages (typically of the order of 10~ 5 V/K) 
at the connectors. Although they contribute only on the sub-ppm level they can add up to a signifi- 
cant uncertainty. Therefore we are using special PTFE LEMO cables (twisted pair with shielding) 
and Cu-Be or Cu-Te connectors while keeping the environmental temperatures stabilised. The 
electric connections are realised with crimping instead of soldering. As for the HV connections 
clean contact and insulator surfaces were maintained. At HV above 10 kV the relative error caused 
by the thermoelectric effect is in the 10~ 10 range and therefore negligible. Instabilities sometimes 
observed at higher voltages of about — 30 kV (used for the measurement of high energy conversion 
lines of 83m Kr, e.g. L 3 -32 line) can be removed by cleaning the HV insulator. However, erroneous 
measurements may still occur due to insufficient electrical connections but should be recognised 
by the redundant monitoring with classical electric techniques and the nuclear standard presented 
here. 

The solid 83 Rb/ 83m Kr sources are intended for continuous monitoring of the energy scale sta- 
bility of the KATRIN spectrometers. This way they should act as an independent stable reference. 
However, in our test measurements with the sources the HV dividers K35 and K65 were used as 
a reference. They were regularly calibrated at the PTB and the drifts of their dividing ratios were 
well known. According to the concept of the continuous monitoring of the energy scale stability in 
KATRIN, the 83 Rb/ 83m Kr source will be biased by — 0.8kV when used at the monitor spectrom- 
eter. This way it will be ensured that the ^"-32 conversion line will match the fixed HV of about 

7 Such a low drift can be safely neglected in the considerations presented later in this work as the observed drifts of 
the energies of the conversion lines were higher by at least one order of magnitude. 




(3.1) 
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— 18.6kV applied simultaneously to the main and monitor spectrometers. However, in our case 
the tested source was grounded and only the HV applied to the spectrometer electrodes was varied. 
The reason for such a configuration was a minimisation of the HV measurement uncertainties. 

3.2 Typical measurement and fit of a conversion line 

The z position relative to the solenoid B centre was common for all the sources. At the beginning 
of the measurement of each electron source the x-y position of the table was found at which the 
maximal amplitude of the conversion line K-32 was observed. Such a procedure ensured that 
each source was individually centred with respect to the spectrometer z axis and that, in turn, the 
electrons emitted from each source experienced the same potential dip across the analysing plane. 
The remaining sources were biased by + 100 V, reducing the electron energy by 100 eV, so that they 
had a negligible effect on the spectrum measured with a given source. In addition, the signal-to- 
background ratio was optimised by tuning the currents in the Earth's magnetic field compensation 
coils which aligned the electron flux inside the spectrometer. 

The actual measurement consisted of scanning a chosen electron line of the 83m Kr conversion 
spectrum. The scan was realised by stepping the HV in a suitable region (typically ±25 V) around 
the voltage corresponding to the energy of the line. The uniform distribution of the measurement 
time in this interval was chosen without further optimisation. The stabilisation time of the HV scale 
due to the initial warmup of the various devices (power supply, divider, digital voltmeter) of the 
order of ~ 12 h was always ensured. Prior to the measurement in each point the HV power supply 
was set to a given step value £/ st ep (requested by the control program) and it was let to stabilise 
for about 5 s. A typical step size was 0.5 V. In addition, after large HV steps, e.g. after switching 
from — 17.8kV (scan of the K-32 line) to — 30.5kV (L3-32 line), a new measurement was started 
only about 30 minutes after the previous scan had finished. Therefore, the complete HV scale was 
always well stabilised in the ppm precision range. 

For each measurement point at a given £/ ste p the measurement real time t r , mean value U meas 
from altogether N low voltage values measured by the voltmeter and a raw ADC spectrum were 
obtained. Typical ADC spectra are depicted in figure U The counts in the regions covering the 
electron peak and the pulser peak were summed up, resulting in peak areas N e and N p , respectively. 
Due to the periodic nature of the pulser the number of events in this peak does not followed a 
Poisson distribution but a binomial one. Therefore, in accordance with [ 5311 the variance reads 



where / denotes the frequency (nominal value of 100 Hz) of the aforementioned pulse generator. 
For the details of the dead time correction see [BOl . 

Each HV region was sweep-scanned in both "down" (absolute value of negative HV decreas- 
ing) and "up" (absolute value of negative HV increasing) directions, while both data sets were 
considered separately as individual integral spectra. One measurement consisted usually of several 
sweeps. Typically, one sweep took about 0.5-4 h, depending on the preset time and the size and 
number of HV steps. As the radioactive decay of the solid 83 Rb/ 83m Kr sources is governed by the 
half-life of 83 Rb, the time duration of one sweep was negligible with respect to r 1 / 2 ( 83 Rb) and thus 
no half-life correction of the counts within the sweeps was necessary. 




(3.2) 
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Figure 6. Typical raw ADC spectra obtained during the measurement of the line of the source Pt-30-2. 
The spectrum marked red was recorded in the condition of "zero" spectrometer transmission. In this situ- 
ation the HV applied to the spectrometer electrodes is rejecting all electrons of the Lj,-32 line and only the 
background — caused by 83m Kr conversion lines of higher energies and by secondary electrons emitted from 
the electrodes — is recorded. The spectrum marked blue was recorded in full transmission, i.e. when all the 
L3-32 electrons are passing the filter HV. The perpendicular lines frame the summation windows in which 
the electron and pulser (around the channel 1700) peaks are summed in each spectrum. The comparison 
of the electron peak areas in the depicted spectra results in the signal-to-background ratio of ~ 18. As a 
consequence of the dead time effects due to the high count rate in the case of the full transmission spectrum 
the pulser peak area is smaller by about 13 % in this case. In the full transition spectrum a double electron 
peak is visible around the channel 1370, stemming from the pile-up effect. 



A typical integral spectrum "dead time corrected count rate in the signal window versus the 
spectrometer retarding energy" is shown in figure [7[ In general the electron spectrum measured 
with the MAC-E filter can be expressed as 

g(qU)=R(E stan ,qU)<g>f(E sta «) +B= (T ®/ trans . loss ) ® (/ el . loss <g> / shake ®V) +B, (3.3) 

where f(E start ) represents the energy distribution of the electrons emerging from the source and R 
denotes the response of the complete system to monoenergetic electrons. The typical Lorentz dis- 
tribution L describing an excited atomic state with the natural width Y is smeared by the Gaussian 
distribution G (broadening due to the Doppler effect and HV fluctuations and possibly due to un- 
resolved chemical shifts), resulting in the Voigt function V = G®L. The term /shake stands for the 
shake-up/off phenomena representing the intrinsic mechanisms of the energy losses of electrons 8 . 
Further, the loss function f el i oss describes the effect of the electron inelastic scattering within the 
source. The inhomogeneities of the electric and magnetic fields 9 can be expressed analytically and 

8 The creation of an inner shell vacancy in the atom — here due to the conversion process — may lead to an additional 
simultaneous excitation or ionisation of outer shell electrons. In the shake-up process a bound electron is excited into 
an unoccupied bound orbitals. During the shake-off process a bound electron can even be ejected into the continuous 
spectrum of free electron states. The original electron, ejected from the inner shell (and thus causing the initial inner 
shell vacancy), is affected by the shake-up/off processes and loses a portion of its energy. 

9 In an ideal MAC-E filter both the electrostatic retardation potential U and the minimal magnetic field strength 6 m ; n 
are invariable along the full extent of the analysing plane. In reality, however, inhomogeneities both in the electric 
potential and magnetic field will be introduced due to the effects of the finite size of the electrodes and the auxiliary coils 
as well as local distortions of the field configuration. 
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Figure 7. Integral spectrum of the K-32 line of the source Au-30 and the fit on the range of [qU\,qU2\ = 
[17825, 17845] eV. (a) The recorded spectrum corrected for the dead time is shown in red. For the sake 
of clarity the error bars of the count rate are multiplied by the factor of 10. The blue line denotes the 
least-squares fit of the convolution of the spectrometer transmission function with a single Voigt pro- 
file. The vertical line at 17 825eV marks the lower limit of the fit range. The Lorentz width was kept 
fixed to r = 2.7 eV, whereas the Gaussian width a was a free parameter. In this case the fit delivered 
a = 0.77(2) eV. The green line denotes the Voigt profile resulting from the fit (including the constant back- 
ground B = 1313(1) counts/s). The peak centroid was fitted as E' — 17 827.827(12) eV. (b) The normalised 
residuals of the fit (A^of = 98, Xi = 1-10) shown in (a) exhibit no significant structure. 



included in the transmission function T, while /trans, loss represents the transmission losses which 
occur typically at high surplus energies as the weak magnetic field in the analysing plane region no 
longer guides the electrons H5H H - Finally, B denotes the recorded background rate. 

For the purpose of the long-term energy stability measurements of the 83m Kr conversion elec- 
trons only the zero-energy-loss peak was fitted on the energy range which excluded the shake-up/off 
peaks, the inelastic scattering of the electrons as well as the effect of spectrometer transmission 
losses. The fit range was selected also on the basis of Monte Carlo calculations of the elastic and 
inelastic scattering processes of conversion electrons emitted from the implanted sources [BOl . In 
this case the spectrum description simplifies to 

,?zero-en.-loss (qU) = T®V +B. (3.4) 
The function g Z ero-en.-ioss was least-squares fitted on the dead time-corrected spectrum with 



-17- 



the help of the MINUIT fit routine H55H built in the ROOT environment H5H1 . For one elastic peak 
altogether four parameters were used: peak centroid E', line amplitude A, Lorentz width T and 
Gaussian width a. Usually, T was set fixed to the literature value for a given conversion line, while 
a was a free parameter. 

3.3 Cross-correlation fit method 

The precise knowledge of the line shape is crucial for a good fit of the conversion line. However, 
a certain statement about the long-term stability of the electron line energy can be made even 
without a detailed knowledge of the line shape (as well as the transmission function shape). The 
method of cross-correlation is often used in signal processing for determining the similarity of two 
signals. Accordingly, two integral conversion electron spectra were compared with each other: 
One spectrum, usually measured with high statistics and a fine energy step, was considered as 
a reference which was then matched to other spectra. The method was utilised in the following 
manner: Firstly, the reference spectrum, in the case of the K-32 line measured with the step of 
0.2 eV, was considered on a range covering the electron elastic peak and background on the high 
energy side. Secondly, such a template spectrum was least-squares fitted onto the new spectrum 
(measured with a usual step of 0.5 eV), while the amplitude A, the spectrum position shift AE' and 
the constant (energy-independent) background B were free parameters. The linear interpolation of 
the template was necessary in order to obtain its y values in the x points (retarding energies) of the 
new spectrum. Finally, the goodness of the fit was verified with %v an d normalised residuals. 

This way a statement about the electron line shift can be easily done as only three free pa- 
rameters are fitted. However, such a method works flawlessly only if the shape of the conversion 
line does not change with time. On the contrary to the vacuum-evaporated 83 Rb/ 83m Kr sources 
where the 83m Kr conversion lines could be easily described [BOl with a single Voigt profile, in the 
case of the implanted sources it turned out to be very difficult to describe the electron spectra rea- 
sonably well with a single peak on the full range. However, omitting the low energy side of the 
peak yielded a good fit, cp. figure [7]. Structures in residuals indicated the presence of a certain 
asymmetry on the low energy side of the elastic peak. It was possible to fit such an asymmetry with 
a second Voigt profile with the same widths Y and a as the main peak. The possible origins of the 
doublet line structure are discussed elsewhere fBOl . Due to certain difficulties of interpreting the 
low energy side asymmetry and its influence on the fitted centroid of the main peak in the classical 
many -parameters fit method the cross-correlation fit method was used instead in order to determine 
the long-term drifts of the conversion lines. 

The cross-correlation method was compared with the classical many-parameters fit mentioned 
above. For the comparison the data of the K-32 line (w4 x 10 3 counts/s) of the source Au-30 
were taken, measured in well defined experimental conditions spanning 26 days. A slight discrep- 
ancy of ~ 4 x 10~ 7 month -1 between the ^-32 line energy drifts obtained by both methods was 
found. This can be ascribed to the fact that slight changes of the conversion line shape in time 
like Gaussian broadening, signal-to-background ratio and doublet structure are very likely to occur 
during the long-term measurements spanning weeks or months. We take this observed difference 
of 0.4 ppm/month as a systematic uncertainty into account when analysing the data with the cross- 
correlation method. 
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4. Experimental results 



4.1 Drifts of 8Jm Kr conversion lines 

At the beginning of the measurement all four implanted sources were placed into vacuum simul- 
taneously. About one month prior to the start of the measurements a thorough bake-out of the 
complete vacuum set-up was accomplished. Moreover, after inserting the sources into vacuum the 
source section was moderately heated (32 h at about 80°C) 10 . The vacuum in the MAC-E filter tank 
was stable in the 10~ 10 mbar range. The actual measurement consisted of frequent scanning of the 
most important 83m Kr conversion lines K-?>2 (kinetic energy of 17.8 keV), Li-9.4 (7.5 keV), 
(30.5 keV) and N 2 /3-32 (32. Ike V) of each implanted source. The K-32 line was measured most 
often, however, the other lines turned out to be very important for additional systematic measure- 
ments. The drifts of the K-32 line energy of each source are depicted in figure The line position 
shifts were analysed with the help of the cross-correlation method with a reference spectrum for the 
data of each source. From the beginning of the measurement the ^-32 line positions of each im- 
planted source exhibited its own unique drift in the ppm/month range. As each source was actually 
unique in its properties the differences in the drifts were not surprising. The individual drifts of the 
sources, corrected for the drifts of the HV divider K65 and the digital voltmeter, are summarised 
in table @. The drifts of the sources Pt-30 and Pt-15 are compatible with zero within la and 2a, 
respectively. The positive drift of the source Pt-30-2 is not compatible with zero which makes this 
source different from its predecessor Pt-30. Considering the relative drifts, the three sources based 
on the implantation into Pt foils exhibited the drifts in the sub-ppm range, while the relatively high 
drift of the source Au-30 amounted to 2.39(19) ppm/month. This observation could be attributed 
to the difference between the noble metals of gold and platinum. The drifts of the other 83m Kr lines 
were studied as well — much less frequently than the K-32 line, however, and thus with a lower 
statistics — and the absolute drift values were found compatible with the drifts of the ^-32 of the 
individual sources. 

Figure |9| shows a longer period of data taking with the four sources. The vertical dashed line 
in figure |9| denotes the vacuum breakdown which occurred during the measurement (i.e. all the 
vacuum sections were affected) on day 26 as a result of the power outage in the whole institute 
building u . The vacuum worsened from 10~ 10 up to 10~ 4 mbar. The sudden drop of the K-32 
line positions of all sources is clearly visible in the plot. The weighted average of the K-32 line 
shifts of all the sources amounted to —0.149(22) eV. In the case of the source Pt-30-2 it was pos- 
sible to determine the sudden drop of all four most important 83m Kr lines (Li-9.4 to N 2 /3-32) as 
—0.157(42) eV (weighted average), independent of the line energy. These observations can be ex- 
plained as an abrupt increase of the work function of the spectrometer electrodes spec (common 
for all the lines of all the sources) upon the residual gas adsorption and a slow recovery of the work 
function with time during further measurements after the power outage. 

10 The procedure of moderate heating was intended to reduce the outgassing of the sources and the complete source 
section. Some influence of heat on the sources cannot be excluded, e.g. the distribution of ions in the host generally 
changes at elevated temperatures. However, the common practice of removing the host lattice damage via thermal 
annealing usually requires 2/3 • T me \ t , where r me i t is the melting point temperature of the host. Therefore, we assume 
that 80° C did not significantly affect the distribution of 83 Rb ions in the Pt and Au foils. 

"Unfortunately, no safety system using some kind of uninterruptible power source was in operation in that moment. 
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Figure 8. The relative positions of the K-32 line of the sources Pt-30 (red), Pt-30-2 (green), Pt-15 (magenta) 
and Au-30 (blue) as a function of time. The linear fit to positions is also indicated, where the Pt-30 data 
were fitted with N dof = 92, X? = 1.27 and similarly Pt-30-2 (JV dof = 311, %t = 1-65), Pt-15 (N dof = 248, 
X? = 1.28) and Au-30 (JV dof = 171, %} = 1.87). The resulting drifts are summarised in table || The x axis 
denotes the relative time in days since 27.07.2009. The y axis shows the line energy shift AE' relative to the 
reference spectrum. The full span of the y axis of 0.2 eV represents the portion of ~ 1 1 ppm considering the 
line energy of 17.8 keV. 




relative time [day] 



Figure 9. The relative positions of the 7^-32 line of the sources Pt-30 (red), Pt-30-2 (green), Pt-15 (magenta) 
and Au-30 (blue) as a function of time. The x and y axes have the same meaning as in figure 0. The linear fit 
to positions on the range of days 0-26 is also indicated. The full span of the y axis is 0.35 eV (~ 20 ppm of 
the K-32 line energy). The vertical dashed line on day 26 denotes the vacuum breakdown. The chain-dotted 
line on day 35 denotes the test of deliberate venting of the source section. 



In order to verify the hypothesis that the sudden shift of the K-32 lines occurred due a change 
of spec several tests of deliberate venting of the source section were carried out. The vertical chain- 
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Table 2. Overview of the relative drifts of the K-32 line of the four implanted 83 Rb/ 83m Kr sources. The 
drifts di were determined from the analysis of the data collected in days 0-26 only (when the experimental 
conditions were stable). The drifts bi stem from the joint analysis of all the K-32 line data collected in days 
0-46 (including the shift of the work function spec due to the vacuum breakdown). 



source 


lin. drift d t (days 0-26) 

abs. rel. 
[meV/month] [ppm/month] 


lin. drift b{ (full data in days 0^16) 
abs. rel. 
[meV/month] [ppm/month] 


Pt-30 


-5(6) 


-0.27(32) 


-17(3) 


-0.96(17) 


Pt-30-2 


11(3) 


0.63(18) 


18(3) 


1.02(16) 


Pt-15 


4(2) 


0.23(14) 


3(3) 


0.19(14) 


Au-30 


43(3) 


2.39(19) 


47(3) 


2.64(16) 



dotted line in figure ^ represents one of such tests. The vacuum in the spectrometer vessel was kept 
intact (closed by the gate valves) while the source section was vented with air up to atmospheric 
pressure for about 5 minutes. Right after the pump-down back to the 10~ 10 mbar range — without 
any further bake-out of the source section — the stability of the conversion lines was examined. 
The observed shifts of the K-32 line position were compatible with zero: 8(9) meV (source Pt- 
30), 7(4) meV (Pt-30-2), -3(5) meV (Pt-15) and -6(12) meV (Au-30). The tests with deliberate 
venting of the source section showed that the implanted sources are reproducibly stable against 
abrupt changes of vacuum conditions in their vicinity. The same conclusion was obtained from 
tests carried out in other time periods (of the investigations of the solid 83 Rb/ 83m Kr sources) which 
were not influenced by any vacuum breakdown. 

We clearly see that after day 26 the line position of all four lines comes back to the original 
position supporting our assumption of a recovery of the work function spec - In order to disentangle 
this recovery of the spectrometer electrodes from the drifts of the implanted sources a following 
approach was chosen: The K-32 data of each source were assumed in the form 

y i (t) = \ ai+bi ' t f ° r?<?0 ' (4.D 

[cii + bi • t + A0 spec ■ exp [ - (t - t )/r] for t > t . 

where t denotes time in days (to = day 26), are the parameters of the linear drift of the ^-32 
line position of the given source i, A0 spec represents the shift of the spectrometer work function, 
common for all the sources, and X stands for the "lifetime" of the exponential recovery of spec - 
The definition of such a function assumes that the drift of spec was equal to zero before the vacuum 
breakdown occurred. In addition, it is assumed that the drift of the given source was not altered 
by the vacuum breakdown and that the positive drift observed in the time period of days 26-46 is 
a combination of the "original" drift of a given source and the recovery of spec - The latter one 
was anticipated to have an exponential form as a long-term effect of stabilisation of experimental 
conditions. 

A joint fit of the full data set depicted in figure was performed: The function from eq. Q4.1D 
was used for each source i and altogether 10 free parameters were fitted. Four pairs of the pa- 
rameters ai,bi corresponded to four implanted sources, while A0 spec and % were common for all 
sources. The fit is illustrated in figure [RJ The common shift A0 spec was fitted as — 146(2) meV, in 
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relative time [day] 

Figure 10. The least-squares fit (A^j f = 1530, Xi — 0.96) of the combination of four functions defined in 
eq. Q4-.ll ) to all data of the K-32 relative line positions of the four implanted sources. In each source data set 
the sudden shift A0 spec , caused by the vacuum breakdown, is visible as a negative "jump". 

a good agreement with the value of —0.149(22) eV determined from the shifts of individual K-32 
lines. The lifetime X was determined as 68(5) d which illustrates the need for a thorough bake-out 
of the complete system immediately after such a vacuum breakdown accident. The linear terms b[ 
[ppm/month] are listed in table |[ The comparison with the linear drifts determined from the fits of 
the individual data sets on the range of days 0-26 shows that both approaches deliver comparable 
results. 

4.2 Systematic effects 

At this point the various systematic effects of the long-term measurements of the conversion elec- 
trons energy stability shall be briefly discussed: 

1 . the stability of the MAC-E filter energy resolution 

Equation ( |Z.2| ) describes the dependence of the energy resolution on the magnetic fields. The 
stability of S max was monitored by measuring the currents after energising the supercon- 
ducting solenoids and shortly before de-energising them and this way B max was found to be 
stable at the 0.1 % level. The stability of B m [ n (typically 5 m j n ~ 3 G) was estimated by a four 
days measurement outside of the spectrometer vessel to be better than 0.5 %. Both instabili- 
ties can be safely neglected because they do not affect the observed line position at the ppm 
precision. 

2. the stability of the currents in the auxiliary perpendicular coils 

The influence of the currents in the auxiliary perpendicular coils on the energy and shape of 
the conversion line was studied. The parameters of the conversion line (position, amplitude 
and background) showed a complicated dependence on the currents, however, it became 
clear that the auxiliary perpendicular coils affect the line shape in the same way as moving 
the source in the x-y plane. Nevertheless, according to the power supplies specifications 
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an instability of the coil currents has a negligible influence on the position of the K-32 line. 
Finally, any malfunction of one of the coils would be immediately observed as a significantly 
shifted spectrum or as a spectrum with a reduced amplitude-to-background ratio. 

3. the position of the source and detector 

The influence of the source position on the K-32 line was studied with the help of the source 
Pt-30. The reproducibility of the source positioning in the x-y plane was verified. During 
this the K-32 line spectra, recorded after repeated alignment of the source by means of the 
x-y-z table, were compared. The set-up was found to be sensitive to the source misalignment: 
Depending on the direction in the x-y plane, the 1 mm misalignment could cause the shift of 
the K-32 line of 33(3) meV. As the movement precision in the x and y directions amounted 
to 0. 1 mm, the effect of the source misalignment with respect to the drift of the conversion 
line energy could be neglected. The z position of the source with respect to the centre of 
the solenoid B influences the magnetic field B s which, in turn, determines the maximal angle 
(?start of the acceptance cone for the emitted electrons. This z position was fixed during the 
measurement. The position of the detector with respect to the centre of the solenoid A was 
chosen in the following way: The z position was chosen as close to the solenoid A centre as 
possible 12 and the detector was positioned in the x-y plane so that the observed background 
was minimal. Thus, the detector was assumed to lie on the spectrometer z axis. The detector 
position was kept constant during the measurement. 

4. the stability of the spectrometer work function 

In the course of the measurement it became obvious that the stability of spec is a key factor 
which could not be monitored in some way at that time. The only means of controlling 
(minimising) the drift of spec was a thorough bake-out procedure performed at the beginning 
of the measurement. In addition, the simultaneous use (without breaking the vacuum) of the 
four sources allowed to obtain an indication about a possible drift of <p spec as this term is the 
same for all of the sources. 

5. the drift of the dividing ratio of the HV divider 

The drift of M(t) was taken into account in the analysis via the linear term m, cp. eq. fl3.1|). 
The uncertainty of the drifts of the order of O.lppm/month seems to be the limit of the 
systematic uncertainty of the presented method relying on the long-term stability of the HV 
dividers, however, the dividers proved to be more than sufficiently stable for our purpose. 

6. the drift of the scale factor of the digital voltmeter 

On the basis of regular calibrations (once per 1-2 days) the drift of the scale factor K{t) 
was determined as k = 0.39(9) ppm/month. The scatter of the calibration points around the 
least-squares fitted linear function amounted typically to 0.3 ppm. The DC 20 V range of the 
voltmeter was calibrated at — 10 V with the help of the 10 V DC reference. 

7. the stability of the DC 10 V reference 



Due to mechanical obstacles it was not possible to place the detector in the magnetic fields higher than 1.06 T. 



The stability of the DC 10 V reference is specified by the producer as 3ppm/year [231, 
however, the drift of the device used in this work was determined as < 0.3ppm/year \UTW . 
Such a low drift can be safely neglected in our case. 

8. the stability of the scale of the ADC module 

Any unrecognised shift of the ADC scale would directly influence the number of recorded 
electron events which is then used for the construction of the integral electron spectrum. 
However, the summation windows in the ADC raw spectra were always chosen wide enough 
(see figure ^) so that a small shift of the ADC scale would not matter. 

9. the stability of the pulser frequency 

The frequency of the research pulser was regularly (weekly) checked and a trend of stabilisa- 
tion (from ~ 102.5 Hz up to ~ 105.0Hz) was observed during the time period of 7 months. 
In the analysis the mean value was taken into account. The effect of the pulser frequency in- 
stability translates into the change of the fitted amplitude and background of a given integral 
spectrum but it does not effect the fitted line position. 

An attempt was made to critically assess all known systematic effects which could possibly 
affect the presented long-term measurements. Considering the aforementioned effects it seemed 
appropriate to estimate safely the limit of the experimental systematic uncertainty of the line 
K-32 energy drifts as 0.5 ppm/month. Adding this value quadratically with the discrepancy of 
0.4ppm/month stemming from the use of the cross-correlation method the overall systematic un- 
certainty of the drifts reads 0.64 ppm/month. 

5. Discussion and outlook 

With the aim of developing a stable monoenergetic electron calibration source we investigated a 
novel electron source based on 83 Rb / 83m Kr implanted into high purity platinum and gold foils. The 
conversion lines of 83m Kr are suitable for calibration purposes in the KATRIN experiment, namely 
the kinetic energy of the conversion electrons in the K-32 line amounts to 17.8 keV which is by 
only 0.8 keV lower than the tritium /3 spectrum endpoint. At Mainz MAC-E filter we tested the 
stability of four implanted 83 Rb/ 83m Kr electron sources and we showed that they reach the ppm 
range. The long-term stability of the kinetic energy of the conversion electrons was quantified with 
the notion of the relative linear drift [ppm/month] of the given conversion line of the given source. 

The stability of the electrical devices used in this work, namely the high precision HV divider, 
the HV power supply and the digital voltmeter, allowed to extract the electron energy drifts from 
the data. The stability of the spectrometer work function spe c turned out to be one of the key factors 
affecting the measurements. Unfortunately, the set-up used in this work did not allow to monitor 
the stability of spec in an independent way. The only means of controlling (minimising) the drift 
of S pec was a thorough bake-out procedure performed at the beginning of the measurement. In the 
time period of 26 days the drift of spec can be safely neglected and the observed drifts of the four 
implanted sources (see table [I]) can be ascribed to the solid state effects in the sources affecting the 
electron binding energy. The systematic drift uncertainty of 0.64 ppm/month was determined on 
the basis of safe estimates of various effects. All the three Pt-based sources fulfilled the stringent 
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demand of ±1.6ppm/monfh on the energy stability 8E/E set by the KATRIN project. The source 
Au-30 exhibited non-zero drift and one can speculate about the reason for such a relatively high 
drift: Due to the well focused ion beam the implanted dose of 83 Rb might be too high and serious 
damage of the lattice could have occurred. However, the same holds for the source Pt-30-2 where 
the drift was low. This may reflect the differences between the elements of gold and platinum. 

The implanted sources were found to be superior to the vacuum-evaporated ones (which are 
reported in EBD) concerning the resistibility of the conversion line energies to vacuum conditions. 
Several vacuum breakdowns occurred during the measurement campaign which caused severe 
temporal worsening of the vacuum from 10~ 10 mbar up to 10~ 4 mbar. In general, these vacuum 
breakdowns resulted in negative shifts of the conversion line energies measured with the implanted 
sources. The tests of deliberate venting of the source section separated from the spectrometer ves- 
sel made it possible to state that the shift of spec is responsible for the observed shift of conversion 
lines energies. It was found that spec was recovering exponentially after the vacuum breakdown. 
The implanted sources exhibited no shift of the conversion line energy even after exposure to air. 
This feature was reproducibly tested three times with the precision of 4 meV. Any vacuum break- 
down should certainly be avoided in the future measurements — this holds not only for the monitor 
spectrometer, but obviously for the whole KATRIN set-up. The fluctuations of the work function 
will be monitored with the help of dedicated methods on places in the KATRIN set-up where the 
stability of the work function is important regarding the observable m (v e ). 

In addition, the proper HV connection between the main and monitor spectrometers is of a 
great importance with respect to the stability of the energy scale in KATRIN: Any unrecognised 
thermoelectric or cable effect could cause a systematic shift in the monitoring scheme, similar 
to the shift caused by the work function change. Therefore, special care has to be taken of all 
cables, feedthroughs and connectors on the HV as well as low voltage side. Recently, different 
cabling schemes and various locations of the KATRIN HV dividers and HV power supply were 
investigated with respect to the stability of the 83m Kr conversion lines measured at the monitor 
spectrometer at the KIT. No significant change of the stability was observed even when a HV cable 
of length of several tens of meters was used [157I1 . Still, the influence of ground loops and RF ripple 
has to be investigated in the near future when the main spectrometer is running and connected to 
the monitor spectrometer. 

It is worth mentioning that the source Pt-30 exhibited the aforementioned drift, compatible 
with a zero drift, after almost 9 months since its production at the ISOLDE facility. In this sense 
such a solid 83 Rb/ 83m Kr source is indeed very well suited for the long-term application in the 
KATRIN experiment. The four implanted sources showed different drifts probably due to the fact 
that the sources were produced with different parameters. Therefore no strong statement can be 
done regarding the reproducibility. This will be the subject of further studies. Still, it is promising 
that three different sources based on implantation into platinum exhibited the stability compatible 
with the requirements of KATRIN. 

The electron count rate recorded in the K-32 line amounted up to « 10 4 counts/s in the case 
of the source Pt-30-2. Therefore, the concept of dead time correction utilised in this work was an 
important part of the data treatment procedure. The amplitudes of the conversion lines generally 
followed the radioactive decay of 83 Rb with the half-life of 86.2 d, however, significant deviations 
from this trend were observed. The data analysis revealed that the fitted half-life was actually lower 
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(7\ /2 (Pt-30) = 77.5(4) d, 7\ /2 (Pt-30-2) = 78.8(3) d, r 1/2 (Pt-15) = 77.2(3) d and r 1/2 (Au-30) = 
84.3(3) d) than the literature value. However, here it should be noted only the electrons emitted 
from the sources without any energy loss were taken into account this way. The lower than expected 
half-life can point towards the fact that 83 Rb diffusion processes take place in the sources. In 
addition, it cannot be excluded that 83m Kr retention of the sources changes in time. 

Although very promising results were obtained already with the first samples of the implanted 
sources, the following details of the production technique could be improved. Stated are also 
suggestions concerning the future sources which could clarify the open issues. The implanted dose 
Q w 10 14 ions/cm 2 of the current samples shall be lowered by at least one order of magnitude in 
future. This way the lattice damage induced by the 83 Rb ions should be reduced and possibly the 
line shape shall be less asymmetrical. In order to achieve the 83 Rb activity sufficient for long-term 
application of the source, the option of sweeping the ion beam over the substrate can be utilised at 
the ISOLDE facility. Thermal annealing of the sources after the 83 Rb collection shall be tested (it 
is also a common procedure to thermally anneal the substrate even prior the implantation). It was 
realised that according to simulations about 5-10% of 83 Rb ions do not reach the metallic state 
during the implantation, but they end up in the rest gas layer covering the metal foil. Although 
this fact was not actually proved to be important in the presented measurements, the surface of 
the metal foil could, in principle, be cleaned by simple means of ion sputtering or plasma etching. 
Such a procedure would ensure that virtually no 83 Rb atoms are present in the non-metallic state. 
Very promising results were obtained also with the source Pt-15 which, furthermore, profited from 
reduced portion of inelastically scattered electrons. The option of implantation at low energies (« 
5-10 keV) should be tested as it would increase the number of zero-energy-loss electrons, useful for 
the energy scale monitoring. As the substrate the polycrystalline platinum foils of high purity seem 
to be sufficient for our purpose. One could also attempt to implant 83 Rb into the HOPG substrate 
for further reduction of the inelastic scattering of electrons. The use of gold foils does not seem 
promising with respect to the non-zero shift of the K-32 line energy. Finally, the presence of the 
doublet line structure in the conversion electron spectra remains an issue for the precise description 
of the line shape [fJQl . Further measurements with thermally annealed implanted sources may 
clarify the origin of the doublet line structure. 

In summary, the 83 Rb/ 83m Kr sources based on implantation into platinum foils provide elec- 
tron calibration sources with a very high stability in time of the order of 0. 1 ppm/month over several 
months. The achievable count rates are high enough to check the stability of the energy scale with 
15 meV precision within several minutes. Therefore such a source is ideally suited as calibration 
source in the monitor spectrometer of the KATRIN experiment. For the sake of redundancy at least 
two implanted 83 Rb/ 83m Kr sources should be placed in the source section simultaneously and their 
conversion lines should be measured on a regular basis. We also highly recommend to measure 
regularly other conversion lines besides the K-32 line. We will report about the investigation of 
the full 83m Kr conversion electron spectrum at another place [5811 . This procedure proved to be 
very useful for checking the stability of the HV scale and the effects of abrupt changes of <j) spec . 
However, since the source contains a parent isotope 83 Rb with the half-life of 86.2 d an application 
in the KATRIN main spectrometer is not foreseen in order to avoid any risk of contamination in 
case of failure, e.g. by an electric discharge. Its compactness allows the implanted source to be 
also used in other applications, e.g. for the tests of alignment and transmission properties of the 
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individual components of the KATRIN main beam line. 

The presented implanted electron sources can be of interest in other applications where a 
source of electrons of precise and ultra-stable energy is necessary. The energy resolution should 
amount to AE < 5 eV in such applications in order to separate the zero-energy-loss electrons from 
the electrons undergoing the inelastic scattering as the latter ones can be influenced by the changes 
in the source surface. 
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